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Abstract 

This article describes a general method for computing automorphic 
forms using Voronoi-type summation formulas. It gives a numerical 
example where the technique is successful in quickly finding a cusp 
form on G-L(3, Z)\GL(3, M), albeit one whose existence was already 
known as a Langlands lift. 

MSG 2000 classification codes: 11F55, 11M41, 11Y99 

The purpose of this note is to describe a method for computing general 
automorphic forms. I have carried out only limited computational tests so far, 
and have not discovered any new automorphic forms using it. However, the 
method does identify some lifted cusp forms on GL{3), and until recently was 
the only general method to compute an automorphic form on a higher rank 
group. It generalizes the methods of Hejhal (e.g. [8]) and others for Maass 
forms on GL{2), but does not require Whittaker or Bessel functions. Aside 
from trying to compute some of the same objects, this method is otherwise 
unrelated to the cohomological methods developed by Ash and others which 
use geometric data to compute special types of automorphic forms. 

I shall begin with a highbrow version of the method, and later explain its 
concrete manifestations and how they relate to existing methods. Every cus- 
pidal automorphic representation for T\G, where G is a arbitrary reductive 
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Lie group and T an arithmetically defined subgroup, has associated automor- 
phic distributions as in [T^. These are viewed as objects in C~°°{N), the 
distributions on the maximal unipotent subgroup N of G. Furthermore, they 
satisfy invariance properties, in particular under A^r = F fl A^. Fourier ex- 
pansions for L^[Nr\N) and C~'^{Nr\N) can be computed in an algorithmic 
way using tools from representation theory, for example using [1] or Kirilov's 
method of co-adjoint nilpotent orbits. The invariance under the rest of the 
group F then provides many nontrivial (distributional) identities between 
these Fourier components. So far this has been carried out for GL{3) and 
GL(4) (see, for example, [l2l[T3]) and Sp{A) ([6]). The full set of identi- 
ties characterizes automorphy (a kind of converse theorem), and when these 
identities are integrated against test functions, they yield identities involv- 
ing invariants of the automorphic representation (such as Hecke eigenvalues). 
The main idea of this note is to use these relations to numerically identify 
automorphic forms. 

The integrated identities for G = GL{3, M) and F = GL{3, Z) were worked 
out in detail in [12], where they are the analogs of the classical Voronoi 
summation formulas for GL{2). These GL{2) formulas have the form 

^ane{-na/c) f{n) = \c\^y^ e{na/c) F{n/c^) (see [El (1.12)]), 

(1) 

where a/c is a fraction in lowest terms, aa = 1 (mod c), a„ are the coefficients 
of the form, e{u) = e^'^*", / is a test function, and F a type of transform 
of it (for details please see the statement in [12] )• When / is chosen to be 
the appropriate K-Bessel function, the relation ([1]) amounts to the equality 
between the Fourier expansions of a Maass form evaluated at two F-equivalent 
points in the upper half plane. The generalized Voronoi formulas are similar 
equalities, but for specific Fourier components on Nr\N, and more general 
vectors in the automorphic representation (e.g. not just the spherical vector, 
which corresponds to the JT-Bessel function). 

Such a relation, with the i^'-Bessel function, has been the starting point 
for many numerical computations of Maass forms, which seek to solve iden- 
tities such as ([1]) for the values of the coefficients a„. The essential point of 
the method here is that the Voronoi formula still affords this opportunity, 
and in fact with much greater flexibility in that it does not force a single 
choice of test function. Indeed, / can be a certain fractional power times any 
Schwartz function. 
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More importantly, Voronoi formulas have now been established for GL{n) 
in [T2l[I3]- The GL{'i) formula has the following form: 



^aq^ne{-na/c) f{n) = \c\^^Y^S{qa,n;qc/d)F(^'\ 

n^O d\cq n=/=0 ''^ ' \^ 1 / 



(2) 



where S{a,b; c) = J2x€{z/cZ)* e{ "'^~^^^ ) is the usual Kloosterman sum, and F 
represents a different transform of / than in ([1]). The Fourier coefficients on 
GL{3), as they are on other higher rank groups, are indexed by more than 
one integral parameter, though the L-function is determined by the data 

Thus one can simply seek to apply the existing coefficient-solving methods 
to formula ([2]). This seems to be advantageous compared to long-standing 
proposals to explicitly compute GL{n) Whittaker functions: for one thing, 
it is uniform over all representations, and again subsumes the Whittaker 
approach. Thus it cannot be any less successful, though choosing appropriate 
functions may still be a challenge. Also, the Fourier expansion in higher rank 
involves more complicated sums than it does for classical modular or Maass 
forms, and so it is appealing to isolate individual Fourier components of it 
like distributional identities such as ([2]) do. The coefficient-solving methods 
also include techniques to hone in on the correct eigenvalue of form when it 
is not known in advance, which in principal apply here also. I have not been 
able to investigate this in detail, however. 

I have attempted to carry out the method in a simple case, of cohomo- 
logical cusp forms on GL(3, Z)\GL(3, M). The reason for doing so is that 
the representation type is pinned down exactly, and the coefficients are also 
expected to have integrahty properties. Indeed, for many purposes cohomo- 
logical forms are the appropriate analog of discrete series/holomorphic forms 
on GL{2). I wish to stress that this is done simply to avoid precision prob- 
lems, and that the technique is in principal general. The simplest case of ([2]) 
occurs when there is no additive twist, corresponding to a = and c = 1: 

^anf{n) = ^a;;F{n), (3) 

where again / and F are a transform pair of functions. This formula can also 
be proven by writing the lefthand side as a contour integral in terms of the 
Mellin transform of /, ^ J.^^ ^^^L{s)Mf{s)ds, and applying the functional 
equation of the L-function L{s) = 'Y^ann~^. Formula ([3]) is thus essentially 
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n 


Approximate Value 


Exact Value 


2 


-0.718812 


-f ^ -0.718750 


3 


-0.641492 


^ -0.641518 


4 


1.235340 


11 ^ 1.23535 


5 
6 


-0.522150 
0.460477 


1019969 ^ n coooo^ 
1953123 ~ U.dZZZZ^ 

III ^ 0.461091 


7 

8 


-0.854928 
-0.418517 


-iHI - -0-858212 
-ill ^ -0.40451 



Table 1: Coefficients of a cohomological cusp form on 6*^(3, Z)\GL(3, R). 
This form is the symmetric square of Ramanujan's A cusp form. The compu- 
tation took 8 minutes on one core of an Intel(R) Core 2 CPU system running 
at 1.5 GHz, with 2 GB of RAM, slower than the methods of 

an equivalent form of the approximate functional equation. One could also 
of course attempt to use the additive twists e(— na/c) present in ([2]) to gain 
further information. In the experiments f{x) was chosen to be a polynomial 
times a Gaussian, and F was computed using [12l (5.27)] - more accurately, 
approximated by a partial sum of its residues since a closed form expression 
for F was impractical. An important difference between higher and lower 
rank groups is that more coefficients are needed as c grows, owing to the 

growth of the twisted conductor of a degree d L-function, such as one 
coming from GL{d). The smaller growth for GL{2) is an important ingredi- 
ent in Stark's Hecke operator method [TB] (see also [HI pp.220-224]), which 
accordingly appears difficult to generalize to higher rank. 

The experiments were carried out to search for cusp forms whose repre- 
sentation type agrees with the symmetric square lift of Ramanujan's A cusp 
form of weight 12. Here one knows via earlier computations (originating 
in studying the group cohomology of SL{3,Z) - see [1], for example) that 
there is only one such form to be found: the global lift of A. Locking in 
the value ai = 1 and using 02, ... , 050 as 49 unknowns, 49 different Voronoi 
equations were created. The truncation of identity ([3]) to its ffist 50 terms 
on each side was then solved for 02, ... , a^o using matrix inversion. The ffist 
few coefficients found this way are listed in Table [H along with their correct 
values (which can be independently calculated using the description of the 
symmetric square lift). One could further use the fact that the correct values 
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of a„n^^ are integers to eliminate precision errors (in this atypical situation). 

It is natural to ask whether or not there are cohomological cusp forms 
for GL{3,Z)\GL{3,M.) aside from these symmetric square lifts from GL{2). 
They do exist for general congruence subgroups, but none are known at full 
level - this despite a serious numerical search by Ash and Pollack [1], who 
consequently conjectured that none exist. (See also [5J, who propose another 
explanation for the paucity of these forms.) 

In the case of Maass-type (i.e. spherical) cusp forms on GL{3), one does 
not know the eigenvalue to start with, and thus must apply procedures similar 
to those developed for the GL{2) case to find it along with the coefficients. 
Though there are again symmetric square lifts of GL(2) Maass forms, it was 
proven in [TT] that there do exist many nonlifted forms: in fact the lifted ones 
comprise only a minuscule portion of the spherical automorphic spectrum on 
GL{n) [13]. Nevertheless, none have been identified to date, and finding 
them remains an important challenge. 

Recent Developments 

Shortly after this note was first distributed on the arXiv, significant progress 
was made in computing Maass-type automorphic forms on GL{3). The meth- 
ods used are similar to those described above in that they use an integrated 
form of the functional equation, as opposed to trying to use the Whittaker ex- 
pansion directly. However, the success of these methods depends very much 
on some important differences. 

First, Andrew Booker and Ce Bian announced the breakthrough calcu- 
lation of new automorphic forms for SL{3,'Z)\SL{3,M.)/ S0{3,M.) [3]. Their 
method uses the approximate functional equation ([3]), but also for the L- 
function twisted by Dirichlet characters. This is equivalent to using the 
additively-twisted Voronoi formula ([2]), though it is not clear yet whether 
grouping into additive or multiplicative twists is more efficient computation- 
ally. A key point in their work is creating sufficient numerical stability for 
their large system of roughly 10,000 equations (related to [2]). The payoff 
is clear: their results go well beyond the experiment above, which used only 
49 equations. The use of twisted functional equations is motivated by the 
GL{3) X GL{1) converse theorem of [9], which shows that these L-functions 
are sufficient to characterize automorphy. The identity ([2]) shares this fea- 
ture, as do the generalizations to C~°°{Ny\G) Fourier expansions mentioned 
earlier. For this reason, it may be possible to use the distributional identities 
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in [0] to compute cusp forms for 5*^(4), e.g. Siegel modular forms. 

Days later, David Farmer, Sally Koutsoliotas, and Stefan Lemurell an- 
nounced another method to numerically compute automorphic forms on 
SL{3, Z)\SL{3, R)/S0{3, M) [7]. They again use the approximate functional 
equation ([3]), with a wider choice of test functions. However, a crucial point 
is that they force certain expected Hecke relations (such as 0203 = and 
02(25 = Clio), thereby obtaining a nonlinear system of equations. They report 
that this has the effect of drastically lowering the number of equations needed 
down to about 20 for the smallest examples. Despite this reduced number of 
equations, their numerical accuracy still exceeds that of Table [H 

Finally, a third method was introduced by Borislav Mezhericher [TU]. 
This method develops a formula which is simultaneously both a consequence 
of ([2]), and yet a generalization in that it reduces to ([2]) when an auxiliary 
parameter is set to zero. It shares some features in common with Hejhal's 
GL{2) approach, and may be useful for employing techniques that are suc- 
cessful there. 

It should be noted that, in computing spherical automorphic forms for 
GL{n,Z), there are no integrality results that can be relied upon: every pa- 
rameter of these forms, whether they be coefficients or eigenvalues, is thought 
to be transcendental. Thus all methods include an approach to identify the 
correct eigenvalue parameters. 
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